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Abstract 

This paper describes topological kinematics associated with the stirring by rods 
of a two-dimensional fluid. The main tool is the Thurston-Nielsen (TN) theory which 
implies that depending on the stirring protocol the essential topological length of 
material lines grows either exponentially or linearly. We give an application to the 
growth of the gradient of a passively advectcd scalar, the Hclmholtz-Kelvin Theorem 
then yields applications to Euler flows. The main theorem shows that there are 
periodic stirring protocols for which generic initial vorticity yields a solution to 
Euler's equations which is not periodic and further, the L°° and L^-norms of the 
gradient of its vorticity grow exponentially in time. 

1 Introduction 

Knots are an essential ingredient of three-dimensional topological fluid dynamics. Their 
presence as flow or fleld Hnes is a marker of a certain level of complexity in the system. 
A fundamental topological principle is that codimension-two is the knotting dimension, so 
circles can be knotted in 3D and two-spheres in 4D. In a two-dimensional flow, points are 
codimension two, but can they be knotted? Clearly not if we consider the static problem, 
but if we consider their motion, it is plausible that if they get sufficiently entangled during 
their evolution they can have implications for the surrounding flow fleld. 

There is a substantial body of mathematics available to understand this situation, most 
prominently, the Thurston-Nielsen theory. This theory has many aspects and ramifications, 
but the part that is most directly applicable to two-dimensional fiuid dynamics concerns the 
growth of material lines. Specifically, there are certain motions of points or rods in a fiuid 
which imply that the length of a class of material lines must grow exponentially as they are 
advected by the fiuid. Further, the theory gives many methods to detect and/or construct 
such motions as well as algorithms for computing the rate of exponential growth. 

In a two-dimesional incompressible fiuid when a material line is being stretched expo- 
nentially, then under the tangent map some vectors are growing and others are shrinking 
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exponentially. Heuristically, this has implications for a (typical) passively advected scalar: 
its gradients must also grow exponentially. For an Euler flow, the Helmholz-Kelvin Theorem 
says that the vorticity is passively advected and thus its gradients are growing exponentially. 
However, to make these arguments precise, we need conditions on the advected scalar to 
make sure it has nontrivial gradients and more importantly, we must eliminate the possi- 
bility of " unfortunate coincidences" where the gradients often line up with a contracting 
direction in the fluid. 

As is usual in Dynamical Systems Theory, these issues are dealt with by considering only 
generic scalars. A set is considered topologically large or generic if it is dense, G5, ie. it is the 
intersection of a countable family of open, dense sets. Thus with the appropriate norm (and 
thus topology) on the set of scalar functions we only consider those in a carefully chosen 
dense, G^-set. For the case of Euler flows, we use the fact that the initial vorticty determines 
the initial velocity field and so the initial vorticities "parameterize" the collection of Euler 
flows. Our results concern the behaviour of typical values of this "parameter" . 

The second component of this paper studies the situation where a time-periodic stirring 
protocol gives rise to a time-periodic fluid motion. This component does not use the TN- 
theory but shares the two themes of growth rate of material lines and generic hypothesis. For 
a time-periodic fluid motion the presence of a time-periodic passively advected scalar implies 
that the flow is integrable and thus, in particular, the growth rate of advected material lines is 
at most linear. Once again, the Helmholtz-Kelvin Theorem yields an immediate application 
to Euler flows. 

For more information on Thurston-Nielsen Theory see |Thu88j and |FLP91j . for its 
dynamical systems applications see |Boy94| , and for its fluid mechanical applications see 
|BASOO] and |TF06j . For many more references (and pictures of fluids being stirred by pseu- 
doAnosov protocols) see the papers by Thiffeault and Stremler in this volume. The theorems 
in this paper appear in |Boy05| in a somewhat different form; see that paper for details of 
the proofs. 

2 Topological kinematics 

2.1 Basic definitions and terminology 

We first formalize the stirring of a planar body of fluid by moving rods. The fluid region 
is a smooth, one-parameter family of smooth, multi-connected, compact, planar domains 
denoted Mj. In this family the outer boundary is held flxed while the inner disks move. We 
shall always assume time-periodicity of the domains with period one, so M^+i = Mt for all 
t. The moving inner regions are called the stirrers, and they are perhaps permuted each 
period. 

Since we will be using a fair amount of dynamical systems terminology and what is 
called a flow in dynamical systems is called a steady flow in fluid mechanics, we shall adopt 
the terminology that a fluid motion is a smooth one-parameter family of diffeomorphisms, 
0j : Mq — )■ Mt, with 0o = '^d- We may view 0f as Lagrangian fluid displacement map: the 
particle at x G Mq at time is at 0f(x) G Mt at time t. However, at this point we are 



2 



making no assumption that the fluid or its velocity field 



satisfies any particular equations. Because the fluid does not penetrate the moving boundary 
(expressed by 0^ : Mq — )■ Mf) the velocity held satisfles the boundary conditions u-n, = Sj-iij, 
with Bi(t) the motion of the i*'^ boundary. The fluid motion is incompressible or area- 
preserving if V • u = 0, or equivalently, the Jacobian of (f)t is identically one, det{D(f)t) = 1. 

We have assumed that the stirring protocol is time-periodic. A special situation of 
importance is when the velocity fleld is also time-periodic with the same period. In this case 
we can study the time evolution using the Poincare or time-one map, (pi : Mq — >■ Mq = Mi, 
which satisfles 0„ = 0", with the superscript indicating repeated composition. 

2.2 One-dimensional growth rates 

A material fine in the fluid is described by a smooth arc or simple closed curve (sec) 7. The 
main focus here is how the material line grows in length as it is passively transported by 
the fluid motion, which means we are analyzing the length of 0f o 7 as a function of t. We 
will use two different ways to measure the length of the curve, one topological and the other 
metric. 

2.2.1 Metric growth 

Let it{l) denote the length of the curve 7 with respect to some smooth, periodic family of 
Riemannian metrics on the Mt. The metric growth rate of 7 is the growth of 

-met/ A . ^ti(ptOl) 
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as a function of t. There are many methods available to quantify growth rate. The situation 
here is rather simple and we will usually just be bounding the growth above and/or below 
by simple functions, namely, by cA* with A > 1 for exponential growth and by ct for linear 
growth. 

2.2.2 Topological growth 

The topological one- dimensional growth rate is designed to depend on just the topology of 
the stirring process and to be independent of various details of the fluid motion. There are 
two main ideas. First, we restrict consideration to so-called essential arcs and sec which 
truly "see the topology" of the evolving regions. Second, at each moment in time we don't 
compute the actual metric length of the evolving curve but rather we compute the shortest 
length amongst curves with the same topology. Imagine an advecting arc as elastic; at each 
time we will let it shrink back to the shortest length while maintaining the endpoints on the 
same boundary circle. Maintaining the curves topology means that in the shortening process 
it is not allowed to pass through the solid stirrers. 

These ideas are more formally defined using homotopies. Two arcs 71, 72 : [0, 1] — )■ X are 
homotopic in X if there exists a continuous map V : [0, 1] x [0, 1] — )■ X with T{l,t) = 71 (t) 
and r(2, t) = 72 (i)- Thus F gives a continuous deformation from one arc to the other. 



An essential arc is required to have its endpoints on the boundary and two essential arcs 
are homotopic if there is is homotopy F between them with the additional property that 
the homotopy keeps the endpoints on the boundary, so r(s, 0) and r(s, 1) are contained in 
the boundary for all s, Note that by continuity of F this means that 71 and 72 have their 
endpoints on the same pair of boundary circles. In addition, for essential arcs 7 it is required 
that 7 is not homotopic to an arc wholly contained in a single boundary circle. So, for 
example, an arc whose endpoints are on different boundary circles is always essential. 

We let denote the circle and so a sec in X is a continuous, injective map ■y : ^ X. 
Two sec 7i, 72 : S"^ — )■ X are homotopic in X if there exists a continuous map T : [0, 1] x S*^ — )■ 
X with r(l,t) = 7i(t) and T{2,t) = 72(t)- An essential sec is one that is neither homotopic 
to a point nor homotopic to a boundary circle. Thus an essential sec encloses at least two 
but not all of the stirrers. All the arcs and sec shown in Fig. l.(b) are essential. 

In both cases we denote the set of curves homotopic to 7, ie. its homotopy class, as [7] 
and the topological length L*°p(7) of a curve will be the least length among curves in its 
homotopy class 

L*°P(7) ■=mm{i{a): (x G [7]}. 
The topological growth rate of the class [7] is the growth of 

Thus to compute the topological growth rate we evolve the curve forward for time t and 
then measure the least length in its homotopy class. An crucial property for essential curves 
7 is that 

Lr\i) > L'rii). (2.1) 

2.3 Topological entropy 

A fundamental result for two-dimensional iterated C°°-diffeomorphisms is that the maximum 
exponential metric growth rate of arcs is equal to the topological entropy ( |NP93] . |New91] ). 
Note the potential confusion with the current terminology, namely, the maximal metric 
one-dimensional growth is equal to the topological entropy. This entropy is so-called since 
it may be defined using growth rates of distances using a topological metric, in contrast 
to the measure-theoretic (or metric) entropy which requires an invariant measure. The 
variational principle says that the topological entropy is the supremum of the measure- 
theoretic entropies over all invariant measures: however, for curves, the topological growth 
is a lower bound for the metric growth. With all this potential confusion we hopefully clarify 
matters by emphasizing that our main concern here is the exponential growth of the lengths 
of one-dimensional curves. 

In analogy with the result for iterated diffeomorphisms, for a general two-dimensional 
fluid motion 0t define its topological entropy as 

log(L?^^*(7)) 

htop{4>t) = supjlimsup : 7 is a smooth arc}. 

t^oo t 



In accord with equation (2.1), if some essential curve 7 has LfP(7) > CA* for some constants 
C > and A > 1, then /itop(0t) > log(A) > 0. 
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2.4 Isotopy and braids 



As noted in Subsection 2.2.2 the topological growth rate of curves depends only on the rough 
topology of the stirrer motion. This motion, in turn, determines what is called the isotopy 
class of the map 0i. More precisely, two homeomorphisms f,g : Mq — )■ Mq are isotopic if 
there is a continuous family of homeomorphisms ht deforming one to the other with = f 
and hi = g. Given an essential arc 7, applying the isotopy between / and g yields a homotopy 
Hf = ht o '-f between / o 7 and g o 'j- Thus we see that isotopic homeomorphisms give the 
same topological growth rate of an essential curve 7. 

The isotopy class of a stirring protocol can be visualized and in fact characterized by the 
space-time trace of the stirrers or their braid. The collection of braids forms a group and the 
algebra of the braid gives a convenient method to compute properties associated with the 
trichotomy of the next subsection. 

Another comment that will be relevant later is that because our stirring protocols are 
periodic, 0i is isotopic to 0„ for all n G N, even when the velocity field is not periodic. 

2.5 The Thurston-Nielsen trichotomy 

For the purposes of this paper, the Thurston-Nielsen theory classifies stirred fluid motions 
and their isotopy classes based upon the rate of their topological one-dimensional growth: is 
it exponential or linear. 

Theorem 2.1 (Thurston-Nielsen Trichotomy) Let Mt be periodic stirring protocol with 
fluid motion (pf Then either 

(1) PseudoAnosov (pA): there exist constants A > 1 (the dilation) and < Ci < C2 such 
that for every essential curve 7, 

CiX' < L*°^(7) < C2A*, 

and thus htop{(pt) > log(A) > 0. 

(2) Finite order (fo): there exists a constant K > such that for every essential curve 7, 

L'r{i)<Kt. 



(3) Reducible case: (roughly stated) Mq splits into (pi-invariant subsurfaces on which 1. 
or 2. holds. 

We shall call a stirring protocol finite order, pseudoAnosov, or reducible according to its 
TN-type. A critical feature of the TN-trichotomy is that in the finite order and pA cases, 
every essential curve has the same basic topological growth rate. The pseudoAnosov case is 
the most useful for applications because this rate is exponential. 

It is important to emphasize that the TN theory only concerns topological growth rates 
and thus it only gives bounds for the more physical metric growth. Many or perhaps most 
of the fluid motions arising from the same stirring protocol {ie. in the same isotopy class) 
could have a greater growth rates. In particular, a finite order stirring protocol of a real 
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fluid will almost certainly have regions of exponential metric growth while its topological 
one-dimensional growth is always linear. 

It is also important to note that the theorem was first proved in the much broader context 
of isotopy classes of surface homeomorphisms. It is usually framed in terms of the existence 
of a special map, the Thurston-Nielsen representative, which is present in each isotopy class. 
When this special representative map is pseudoAnosov, it has a wealth of nice properties 
including a Markov partition which yields a symbolic model of a mixing subshift of finite 
type. This has many implications, for example, the pseudoAnosov map is ergodic and mixing 
with respect to Lebesgue measure and has a dense orbit and its set of periodic orbits is also 
dense. Handel's Isotopy Stability Theorem shows that most of the dynamical properties of 
the pseudoAnosov map are present in any other map in the isotopy class, though they could 
exist in a "small" invariant set ( |Han85j ). 

3 Passive advection of scalars 

In this section we give two results concerning the behaviour of advected scalars under a 
stirring protocol. The results continue to be strictly kinematic. Given a fluid motion 0^, a 
function a : Mt x M — )• M is called a passively advected scalar if it is constant on trajectories, 
Q;((0t(x)) = q;o(x), or equivalently, dat{(f)t{x)) / dt = 0, where we have written a;t(x) for 
a(x, t). In the language of global analysis one says that at is the push forward of ao and 
writes (0t)=^(ao) = at, with {(t>t)*{ao) = "o o {(f>t)~^- 

Now for any function / : Mq — )■ M we may obtain a passively advected scalar simply 
by defining at := {4>t)*{f), and so in this sense only the initial configuration ao matters. 
For example, if ao represents the initial concentration of dye in a fluid, at = (0t)*(ao) is 
the concentration after time t. On the other hand, in a physical fluid there may be scalar 
quantities of importance which at each time t are computed from the velocity field. In 
this case the advected scalar represents a conserved quantity having physical meaning. For 
example, in two dimensions the curl, cut = V x u, is a scalar which is passively advected in 
an Euler flow. 

3.1 Consequences of a pseudoAnosov protocol 

For a fluid motion determined by a pseudoAnosov protocol we know from the TN-trichotomy 



and equation (2.1) that the metric length of essential curves is growing exponentially fast. 
This means that tangent vectors to these curves must be growing in length exponentially 
under the action of the space derivative of the fluid motion D(j)t. This implies that an 
eigenvalue of Dcpt is growing exponentially. For a passively advected scalar at, since at = a^o 
{(f)t)~^, we have Vat = Vao(-D0t)^^. If the fluid motion is incompressible then det(L>(/)f) = 1 
and so {D(f)t)~^ also has an eigenvalue growing exponentially. Thus as long as there is not a 
unfortunate coincidence where Veto stays aligned with the stable eigen-direction of {Dcpt)"^ , 
we have that |Vat| is growing exponentially. 

There are a number of issues involved with making this arguement rigorous; the locations 
in 7 where there is tangential growth will be shifting in time and the unfortunate coincidence 
could actually happen, especially if one is working in a restricted class of physical interest. 
One way around these difficulties is to take a more geometric global viewpoint. 
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level sets of scalar 



level set of scalar 



Figure 1: (a) Proof of Theorem 3.1 



(b) Proof of Theorem 3.2 



Theorem 3.1 Let Mt be a time-periodic stirring protocol of pA type with incompressible 
fluid motion 0f is a passively advected scalar such that its intial state ao is a generic 

C"^ -function, then there are positive constants c,c' so that 



for all t G M, where X > 1 is the dilation of the pseudoAnosov protocol. 

Here are the main ideas in the proof. First, find a C^-open, dense set Q inside the Morse 
functions (functions with nondegenerate critical points) on Mq so that ao & G imphes that 
cto has a band of regular inverse images which are essential arcs or circles. The pA protocol 
forces a stretch in length by A*. This coupled with area preservation and transport of vorticity 
force the level sets of at = {ipt)^{aQ) to bunch up in the transverse direction, which causes 
llVatlloo — ^ oo like A*. See Fig. l.(a) 

3.2 Time-periodic fluid motions 

Now consider the special situation when the velocity field of the fluid motion (pt is time- 
periodic with period one and so (pi is a Poincare map. In addition, we assume the existence 
of a passively advected scalar that is also time-periodic. This again is rather special, but 
if the advected scalar depends on the velocity field {eg its curl), it is automatically time- 
periodic when the velocity field is. In this situation we show that if the scalar field is typical 
in the appropriate sense, then there is at most linear growth of the metric length of material 
lines. 

The main observation required is that given a time-periodic passively transported scalar 
at = at+i and the Poincare map 0i, the initial configuration ao satisfies for every x 



Thus ao is what is called an integral of motion for the map 0i and so as long as ao is 
sufficiently nondegenerate, its existence precludes chaotic dynamics. 



sup |Vat(x)| > cA* and 
xeA/o 




ao(x) = ai(0i(x)) = ao(0i(x)). 
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Theorem 3.2 Let Mt be a time-periodic stirring protocol with time-periodic fluid motion (pt 
and time-periodic transported scalar at- If the initial state of the scalar ao has finitely many 
critical points (for example, is and nonconstant or is C'^-generic), there exists a constant 
K > so that 

Lr\i) < Kt 

for all arcs and sec 7. Thus the metric one- dimensional growth is linear, and the Poincare 
map (pi has zero topological entropy and zero Lyapunov exponents almost everywhere. 

The geometric heart of the proof is this: the level sets of ao must be preserved by (pi 
and for a generic function ao the contours ao{x.) = c are smooths arcs and circle for typical 
c G M. Thus the dynamics of (pi consists of one-dimensional invariant subsets. Injective 
maps of arcs and circles yield very simple dynamics with no dynamical entropy and using 
this in the entire fluid domain we have only linear metric growth of arcs. 

3.3 Combining the two results 

Putting these two kinematic results together yields an observation of relevance to Euler 
flows. If a time-periodic fluid motion has a generic advected scalar field that depends on the 



velocity field and thus is also time-periodic we know from Theorem 3.2 that there is linear 



metric one-dim growth of all curves (material lines). On the other hand, we know from 



Theorem 3.1 that a pseudoAnosov stirring protocol always causes exponential stretching of 
some material lines. These clearly can't coexist and so if a fluid motion is stirred by a pA 
protocol and has a passively advected generic scalar then either the scalar is not dependent 
directly on the velocity field or if it is, then the fluid motion is not time-periodic. 



4 Euler fluid motions 

4.1 Basic definitions and results 

Now we assume that the velocity field u(x, t) of the fluid motion (pt satisfies the incompress- 
ible, constant density (p = 1), Euler equation 

— = -Vpt, div(u) = 0, 

with slip boundary conditions on the moving boundary. In this case (pt is called an Euler 
fluid motion. Kozonoi has shown that for the class of two-dimensional problems with moving 
boundary considered here there are classical solutions ( |Koz85j ). The results in this paper 
are all predicated on the assumption that there is a global strong solution with the regularity 
of the initial data: our goal is to analyze its dynamics. 

Recall that for two-dimensional, divergence-free velocity fields a classical result says that 
the curl coupled with the circulations around the boundary components and the boundary 
conditions u ■ Ui = Bj ■ rij determine the field completely. In a certain sense we will view the 
collection of Euler fluid motions as being "parameterized" by possible initial vorticities and 
all our results concern dynamics associated with typical values of these "parameters" . 

The Helmholtz-Kelvin Theorem (1890's) is perhaps the most important result about two- 
dimensional Euler fluid motions: an incompressible fluid motion is Euler if and only if its 
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vorticity is passively transported and circulations around all smooth simple closed curves C 
are preserved, or 

d I 

— ® u • ctr = 

for each boundary circle Ci. The preservation of boundary circulation is a necessary feature 
in multi-connected domains. 

4.2 An Exponential Growth Theorem 

Now we consider an Euler fluid motion under the influence of a pseudoAnosov stirring pro- 



tocol. Theorem 3.1 used with the Helmholtz-Kelvin Theorem immediately yields: 



Theorem 4.1 Let Mt he a time-periodic stirring protocol of pA type with Euler fluid motion 
(pf If the initial vorticity Uo is a generic C"^ -function, then there are positive constants c,c' 
so that 

sup |Vwt(x)| > cA* and / |Vwt(x)| > c A* (4.1) 

xeAfo JMt 

for allt e M, where X > 1 is the dilation of the p A protocol. Thus || Au(x, t)||oo = ||Va;t||oo — > 
oo and ||Uf||(72 — )■ oo, all like A*. 

There are numerous results in the stability literature concerning growth of \VuJt\ for per- 
turbations of two-dimensional steady Euler fluid motions. Yudovich ( |Yud74] . |YudOO] ) and 
others have shown linear growth and Arnol'd ( |Arn72] ). Friedlander and Vishik (( |FV92] ) 
and others have shown the importance of exponential growth. The same basic mechanism 
is in play here as in those results, namely, the growth of under the influence of the 



eigenvectors of the tangent map D(j)t as discussed at the beginning of Subsection 3.1 



4.3 Speculations on applications to more general Euler flows 

We make a few very speculative remarks on how the exponential growth induced by a 
pseudoAnosov stirring protocols might have implications for general Euler flow. In particular, 
it provides a topological perspective on a version of the Yudovich Hypothesis /Conjecture: for 
generic initial vorticity a two-dimensional Euler fluid motion satisfies the exponential growth 
in equation (|4l| (c/. |MSYn8] ). 



There are two main ingredients in this topological perspective. The first is that the 
exponential growth of material lines caused by the stirrers could just as well be achieved 
by advecting points in the fluid ( |Bow78j . |GTF05j ). The TN-trichotomy as stated above 
requires these "virtual stirrers" or "ghost rods" to be periodic points. However, it is clear 
that an aperiodic "tangle" of points evolving can also force exponential topological growth 
of a class of material lines. The second ingredient would say that for typical initial vorticity 
a two-dimensional Euler fluid motion always has such a collection of fluid trajectories. If 



both these ingredients were known, virtually the same proof as Theorem 4.1 would yield the 
version of the Yudovich Conjecture just given. While both ingredients seem quite reasonable, 
rigorous results seem a long way off. The topological ingredient is probably tractable but 
determining even a part of the dynamical evolution of a typical two-dimensional Euler fluid 
motion requires tools still to be developed. 
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4.4 A linear metric growth theorem 

Because the vorticity cot is passively transported scalar for an Euler fluid motion as in Sub- 



section 3^ it follows that when the fluid motion is time-periodic, the initial conflguration 
Uo is an integral of motion for the Poincare map 0i. Brown and Samelson ( |BS94] ) ob- 
served using a theorem of Moser that if the integral of motion uq is real analytic, then 0i 



can't have a transverse homoclinic intersection. More generally using Theorem 3.2 and the 
Helmholtz-Kelvin Theorem we have: 

Theorem 4.2 Let Mt be a time-periodic stirring protocol with time-periodic Euler fluid mo- 
tion (pf. If the initial vorticity Uq has finitely many critical points (for example, is and 
nonconstant or is C"^ -generic), there exists a constant K > so that 

Lr\i) < Kt 

for all arcs and sec 7. Thus the metric one- dimensional growth is linear, and the Poincare 
map 01 has zero topological entropy and zero Lyapunov exponents almost everywhere. 

The crucial features of the theorem are the time-periodicity of the Euler fluid motion 
and its two-dimensionality. Thus a similar result holds in any two-dimensional fluid region. 
It is also important to note what the theorem does not say. It does not say that the generic 
initial vorticity gives rise to a time-periodic Euler flow, and it does not say that among initial 
vorticities which do give rise to time-periodic Euler fluid motions linear growth is typical. 
Understanding these issues requires a knowledge of the size of the set of time-periodic Euler 
fluid motions which are not steady. This is ill-understood at present. There are reasons to 
believe that, in fact, generic initial vorticity never gives rise to a time-periodic Euler fluid 



motion and Theorem 4.2 could be useful in proving this. 



4.5 A dichotomy for time-periodic 2D Euler flows 



Theorem 4.2 is similar to and inspired by Arnol'd's result ( |Arn65] ) on steady, three dimen- 



sional C"^-Euler flows which we frame as a dichotomy: 

(1) For typical vorticity the preservation of the Bernoulli function p-\- ||u|p/2 implies that 
the flow is integrable and typical orbits are conflned to invariant tori and annuli and 
the flow thus has zero entropy and zero Lyapunov exponents almost everywhere. 

(2) For the atypical case of a Betrami flow where the curl is parallel to the velocity fleld, 
one may have measures with positive entropy, positive Lyapunov exponents, etc. 



For time-periodic, two-dimensional Euler fluid we know from Theorem 4^ that in analogy 
with 1., for typical initial vorticity there is integrability and linear metric one-dimensional 
growth. For an analog of the atypical case 2., perhaps the simplest situation is when the 
vorticity is constant everywhere. In this case one always has time-periodic solutions for every 
time-periodic stirring protocols using classical potential theory: if C is a given constant and 
(Fo, Fi, ... , Fm) € M"*"*"^ with ^ F, = is a vector of circulations, then there exists a unique 
time-periodic Euler fluid motion (M^, Uj) with Ut = C for all time t and (^-) u dr = Fj, for 
i = 0, . . . ,m. The proof is essentially to solve the Poisson equation = C* for the needed 
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stream function at each time t and these stream functions depend on time-periodic data 
and so are time-periodic the velocity fields V"*"^. 

Now if we assume the stirring protocol is pseudoAnosov with dilation A, the resulting 
fluid motion satisfies L™^*(7) > L*°^(7) > cA* for any essential curve 7. This implies that 
the Poincare map 0i has htop{4>i) > log (A) > 0, ergodic invariant measures with positive 

metric entropy and thus positive Lyapunov exponents and more Thus one dichotomy 

for time-periodic solutions to Euler's equation in a time-periodic moving domain is: 

(1) For typical initial vorticity there is integrability and the linear growth theorem above. 

(2) For the atypical case of constant vorticity classical potential theory yields time- 
periodic Euler fluid motions and with pseudoAnosov protocols these have chaotic 
dynamics. 

However, once again we emphasize that the dichotomy is strictly conditional. It says nothing 
about the existence or prevalence of non-steady, time-periodic Euler fluid motions. 

4.6 An energy bound 

Conservation of energy is a fundamental property of Euler fluid motions in stationary do- 
mains. For stirring protocols Mt a simple computation yields that the total energy E — l\\u\\l 
is evolving so that 

dE V- / ^ , 

This reflects the fact that the fluid can do work on the stirrers and vice versa and so the long 
term behaviour of the energy in a stirred Euler fluid is initially unclear especially in light of 
the possible exponential growth of ||Va;t||i. We sketch an argument whose main ideas are 
due to Steve Childress which shows that for time-periodic stirring protocols the energy is 
uniformly bounded as t — )■ 00. 

Assume now that the Euler fluid motion u has smooth initial vorticity cuq and circulation 
around the i*^ boundary equal to Fj. By Helmholtz-Kelvin for all times t, u(x, t) has the 
same circulations and further, \u!t\ < fl :— max(|a;o|). We now assume for simplicity that 
the outer boundary of Mt is the unit circle and so Mt is contained in the unit disk D. 
We further assume that each stirrer is circular with area a. Now fix a time t and suppress 
dependence on it. We decompose u into pieces whose energy will be bounded separately. 

First, extend a; to a function cu on with cu = outside D and cu = Tj/a inside the i*'* 
stirrer. Using Biot-Savart we define 

Thus V has vorticity = a; in M, has circulations Fj, and |v| < O in M where fl — max{0, |Fj|} 
using 

max{— // , ^ , c^yl = 1. (4.3) 

This yields ||vf < ttQ^ 
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Next, let h be the harmonic function on M with dh/dn = — v-n^ on each boundary circle 
Ci and V/i has zero circulation around all the boundaries. Finally, let g be the harmonic 
function on M with dg/dn = Bi on each boundary circle Ci and Vg has zero circulation 
around all the boundaries. 

We now have that v + V/i + Vg has vorticity w, circulations equal to Fj, and the normal 
component of its velocity on the boundary circle Ci is ■ rij and so it is equal to u. Because 
(v + V/i) ■ rij = on every boundary circle Ci and is divergence free we have that (v + V/i) -L 
V/i and (v + V/i) -L Vg working in L^. Thus twice the energy of u is 

||uf = ||v + Wh + V^f = ||v + V/if + llV^f = ||vf - II V/if + llV^f < ||vf + IIV^II^ 

Now we reintroduce the time dependence and do the decomposition of u(x, t) into three 
parts at each time t. By Helmholz-Kelvin the circulations for u(x, t) and thus v(x, t) are the 



same for all t as is the value of Vt. Further, the bound using equation (4.3) is independent of 
t and so for all t, ||v(x, t)||^ < ttVL^. Finally, the data determining g{^.,t) is the position and 
velocity of the boundaries and so g is periodic in t, and we get ||V(7(x, t)|p bounded above 
for all time by its supremum over one cycle, completing the time-independent bound on the 
energy of u(x, t). 
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